A = readmatrix ('NiI2_10VCOOL_spectra.csv')
B = readmatrix('supercontinuum_power_calibration.csv')

m = B(28:99,1)
n = B(28:99,2)*0.1

x = A(1:72,1)
y1 = A(1:72,3)*1e12*0.6./n
y2 = A(1:72,5)*1e12*0.6./n
r = (y2-y1)./(abs(y2)+abs(y1))

hold on
plot(x,y1,LineWidth=1,Color=[0.8 0 0])
plot(x,y2,LineWidth=1,Color=[0 0 0.8])
%plot(x,r,LineWidth=1,Color=[0 0 0])
xlim([500,700])
xticks(500:50:700)
ylim([0,0.11])
yticks([0:0.05:0.1])
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off
box on
set(gcf,'Position',[0, 0, 200, 200])